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Abstract 

The recent experimental and numerical tests of cor- 
rugated nozzles have shown some acoustic and thrust 
benefits relative to traditional round nozzles [1-3]. 

For example, a Bluebell nozzle which was obtained by 
3D nozzle design incorporating a corrugated cross sec- 
tion nozzle shape with a sinusoidal lip line nozzle edge, 
can provide an acoustic benefit up to 4dB with ~17. 
thrust augmentation. In references [1-3], this effect 
was explained as being the result of the corrugated de- 
sign producing more efficient mixing of the exhausted 
jet with ambient air. Based on this argument, the au- 
thors [4] have proposed the application of this concept 
for a centerbody (plug) which can form several vortices 
downstream from the centerbody. 

Several different corrugated designs are proposed 
and described in detail in this paper. The main de- 
sign is a Screwdriver shaped centerbody or plug (SCR) 
which was tested experimentally and numerically. The 
acoustic tests were conducted in the anechoic chamber 
of the Central AeroHydrodynamics Institute (TsAGI, 
Moscow) under Civilian Research and Development Foun- 
dation (CRDF) grant. These experiments have shown 
an essential acoustic benefit of ~10-13*/. with the ap- 
plication of the co- annular nozzles by comparison with 
the reference round nozzle with the same mass flow 
rate. However, the expected acoustic benefits with the 
application of the 4-petal Screwdriver shaped center- 
body were not obtained by comparison with the refer- 
ence axisymmetric centerbody (CON) having the same 
length and the same cross section areas at the same 
distance from the nozzle throat. Moreover, for some 
angles ( 6 = 60° amd 90°) noise increase was observed 
(~l-3'/.)- These tests will be continued with the goal of 
obtaining better acoustic results. In particular, acoustic 
characteristics are hoped t.o be improved by moving the 
centerbody into the nozzle and using penetrable walls 
for the SCR and/or for the main nozzle. Preliminary 
results for such approach are very promising. 

Aeroperformance effects were analyzed numerically. 
The nozzle thrust calculations were based on a full Navier- 
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-Stokes equations solver (NSE), and both full and march- 
ing Euler codes: CFL3D [7], CRAFT [8], and Krayko- 
Godunov [9]. Grid preparation and its optimization 
were conducted using GRIDGEN and our own codes. 
The general conclusion of this numerical analysis is some 
thrust loss with the application of SCR design (~1- 
l .5'/,). But, again, some constructive features of SCR 
design give some promising perspectives for its applica- 
tion in aviation and domestic industries. 

1. A Screwdriver centerbody (plug) geometry 
A Screwdriver body shape was discussed many years 
ago in connection with an application of such a con- 
figurations as the front part of bullets, projectiles and 
rockets. This shape provided lower drag and higher sta- 
bility of forward flight for rotating projectiles. These 
surfaces belong to a family of shapes formed out of by 
rectilinear intervals joining corresponding points of two 
different closed curves in space. A Screwdriver body 
uses a circle as an initial curve and one or several cross- 
ing rectilinear intervals as an end curve. Usually, these 
intervals are symmetrically located relative to a body’s 
axis of symmetry. In subsequent discussion, several dif- 
ferent modifications of Screwdriver shaped centerbody 
surfaces will be described. 

1.1 Rectilinear Screwdriver Shaped Sur- 
faces. We use a cruciate Screwdriver shaped center- 
body. A 4-petal Screwdriver centerbody is shown in 
Figure la. Here a convergent nozzle 1 contains a cylin- 
drical centerbody and downstream, after the throat, 
it transfers to a Screwdriver centerbody. The front 
view from downstream is shown in Figure lb. There 
are several geometric parameters which define the cen- 
terbody: number of petals, petal size and centerbody 
length. The axisymmetric centerbody with conical or 
optimal contour in a meridional plane can be taken as 
a baseline centerbody for comparison and for definition 
of Screwdriver centerbody efficiency. Figure 2 demon- 
strates four different modifications of Screwdriver cen- 
terbodies based on the conical centerbody. These shapes 
correspond to a baseline conical centerbody (a), and 
2(b), 3(c), 4(d) and 6-petal (e) contcrbodies. The sur- 
face sketch, coordinate system, nomenclature and sur- 
face equation for a region between t.wo planes of sym- 
metry are illustrated in Figure 2f. These surfaces are 
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formed as follows. The 45 c, -arc of the cylindrical portion 
and vertical interval x t , z° are divided by I subintervals 
[a,-,a l + 1 ] and [6,, 6, + 1 ] respectively, where i=0,l,2,...I. 
Then, the corresponding points a, and b, are joined by 
rectilinear intervals and these intervals form the needed 
surface. 

Our hypothesis is as follows. For example, in the 
case of a 4-petal centerbody, four symmetrical flow re- 
gions occur between two neighboring petals. Three- 
dimensional flows form downstream with swirls at the 
local symmetry line as shown in Figure la. These four 
anti-symmetric swirls meet each other after this point 
at some angle and form large vorticities downstream. 
Schematically such swirl flows are illustrated in Figure 
la. This enhances mixing of the jet exhaust with am- 
bient air and can result in jet noise reduction. This 3D 
centerbody shape also favors flow stagnation at the lo- 
cal symmetry line which increases pressure on the lat- 
eral surface of the centerbody. For the same reason, 
such a 3D flow produces a thinner boundary layer than 
the usual axisymmetric flow. As a result we can expect 
the minima] thrust loss and in some cases even thrust 
augmentation through a Screwdriver centerbody appli- 
cation by comparison with a baseline conical center- 
body. However, the final conclusion can be drawn only 
after additional theoretical research, numerical simula- 
tion and experimental tests. 

1.2 Longitudinal Curvilinear Screwdriver- 
Shaped Surfaces. The second option may be more 
promising for designing a Screwdriver centerbody. It 
can help to avoid possible formation of separation zones 
and internal shock waves immediately downstream of 
the throat. This approach serves to smooth the tran- 
sonic nozzle wall and the sharp wall regions at the sym- 
metry planes. Only two simple method of such smooth- 
ness will be shown. In the case shown in Figure 2f, 
longitudinal rectilinear lines are replaced by curvilinear 
lines so that they are a smooth continuation of corre- 
sponding straight lines on a cylindric portion having a 
horizontal tangent at the end of a centerbody, i.e. at 
the rectilinear ending interval x t z° in Figure 2f. Thus, 
these curves should have points of inflection so that thev 
may be given by the two power functions with conjuga- 
tion at these points of inflection. In the simplest case, 
these functions can be written as follows. Let an ini- 
tial point a, and ending point 6, have Cartesian coor- 
dinates x ° 0 = 0, y° c = r 0 cos(<f>i ), z ° 0 = r 0 sm( fc), where 
r c is the cylindrical portion radius and is the polar 
coordinate of the initial point. For uniform splitting 
4 > = (i/I)n/ 4 , z e = r?(l -i/ 1 ). The curve joining these 
points can be described in Cartesian variables as: 

-(f) = -o-a.-f'' 1 , j/(f ) = yo+a v l ri if 0 < £ < £ r 

( 1 . 1 ) 


^(0 - ze+MWr, y(f) = Mi-f)’* »7fr<f<i 

, ( 1 . 2 ) 
wnere Pi,?i.p 2l ?2 are fixed even powers, and is a 
fixed conjunction coordinate of the two power functions 
which can be varied. The coefficients a, , b. ,a y .b y are 
defined from the conjunction conditions: equality of the 
function and its first derivatives so that, for example, 
for Pi = p 2 =p and ?! = q 2 = q these coefficient are: 

- A _ Pf(tc) 


/(fe) = 


(Zp Zg ) 

-p(l - £c)‘ 


(1.3) 


Several different views of such designed sonic noz- 
zles with the 4-petal Scredriver centerbody are shown 
in Figure 3a-e. The equivalent (with equal cross section 
area) axisymmetric centerbody is in Figure 3f. These 
nozzles were designed and drawn at NASA LaRC and 
then manufactured and tested at TsAGI, Moscow in 
Russia. The results will be discussed below in the 2nd 
section . 


1*3 Bluebell Screwdriver Shaped Surfaces. 

The third possible configuration of a Screwdriver type 
centerbody is based on the Bluebell nozzle concept. Let 
the cross section of a Bluebell shaped design decrease, 
instead of increase as in a divergent conical design. In 
this case, the cross sectional area of the baseline round 
conical centerbody decreases from the maximal value at 
the throat to zero at the axis of symmetry. As before, 
we require equality of both areas but increase the corru- 
gation amplitude, for example, linearly from the throat 
to the end. Then we will form a smooth surface with 
cavities and convexities” like a Screwdriver center- 
body. For example, the cross sectional contour of the 
centerbody may be produced by using the same power 
functions as in the previous case for longitudinal lines 
(1. 1-1.3). The normalized variable £ may be replaced 
by a normalized variable rj = <f>/<f> mas , where <f> is an 
angle cylindric coordinate and 4 > max is a half-period 
of symmetry of a full cross section contour function, 
which for a 4-petal Screwdriver is equal tt/4. Thus, the 
dependence of the radial cylindric coordinate r on the 
normalized polar coordinate 77 is: 


r(r)) = r+ - arf , if 0 < 77 < Vc (1.4) 

r(«7) = r_+ 6(1-17)*, if Vc < r, < 1. (1.5) 

Here again, t) c is a conjunction point; r + (x) and r_(x) 
are maximal and minimal radial coordinates in the cross 
section x=const; p and ? are fixed power parameters; 
and the coefficients a and b are defined from conjunction 
conditions, so that: 
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0 = <m( nr )/»/;- 1 , h = pr/( Vr )/( I - , )f )*- 1 . 



g(i)c) = (r + - r_)l[qr) c + p(l - *7c)]- 

In the particular case when a middle point of an interval 
0 < d> < 4>max - n/n* coincides with a point of inflec- 
tion, it is simpler to use the Cosine function instead of 
equations (1.4, 1.5): 

r (^) = r+ + r+ ~ — cosing), 0 <<f>< r/nt 

Further Screwdriver shaped centrbodies may be de- 
signed with more complex shapes which would promote 
the production of vorticity, thereby enhancing mixing. 
For example, if a straight line edge is replace by an arc 
shape interval, designs such as those seen in figure 4a, b 
may be produced. The design in figure 4a incorporates 
petals which are located symmetrically about the lon- 
gitudinal axis. It would be expected that this geometry 
would not lend a great deal to the production of vortic- 
ity. Figure 4b shows a design constructed with adjacent 
petals located anti-symmetrically about the longitudi- 
nal axis. This would be expected to produce angular 
momentum, which may be useful in some application 
in which several nozzles are involved. Rotation of the 
centerbody may be included to further enhance mixing, 
reduce noise, and augment thrust for this and other 

Screwdriver designs. . 

1.4 Chisel Shaped Centerbody. A Chisel 
shaped centerbody is a modification of a Bluebell shaped 
centerbody and its cross section contour is described by 
a periodic piecewise smooth function. Detailed consid- 
eration of such functions applicable for Chisel nozzles 
are in reference [1]. Two examples of such centerbod- 
ies and corresponding cross section contours are shown 
in Figures 5a-d for 8 and 4-petal designs respectively. 
The general idea of this design is to favor a gradual for- 
mation of a swirl flow downstream with minimal thrust 
loss, and possibly in some cases, with thrust augmenta- 
tion. Wing aerodynamic theory helps to construct such 
a three-dimensional centerbody. A gas flowing along 
a vertical wall produces a larger thrust if corrugation 
convexities are contracted to the centerbody end (as in 
the case of an aerofoil or a diamond-shaped wing). For 
the simplest design in Figure 5, dependence of radius on 
the azimuthal angle in the cross section is described by 
a periodic function r = r(<£>) with a period T=2ir/n c . 
into the first period r = r + = const for 0< <p < <Pi and 
<P 2 < <£ < T, and r = r -= consl for < P ^ ^ 2 ’ 
where y?i = 0.5(T - Ay?) and = 0.5(T + Ay?). We 
call a corrugated surface part a “cavity” or a “convex- 
ity” relative to an internal normal to the nozzle wall . 
A cavity depth (or a convexity height) defined by the 
equality Ar — — r_ , increases along the nozzle cen- 

terline from zero at the initial cross section, x = r a , 
to the maximum value at the exit Ar= A, r£. where 
£ = (r_.r,.)/(x, - i,,), i.e. this coincides with the defi- 


nition of a corrugation amplitude coefficient f>. A cav- 
ity (convexity) width, Ay?= <P 2 - fu also linearly in- 
creases (decreases) downstream from zero (maximum) 
at the initial section to the maximum (zero) at the noz- 
zle exit, i.e. 2y?! = T - A<f-T£. If nonlinear func- 
tions r+(*),r_(z) equal zero at the ends of the interval 
0 < x 0 < x e = 0, the centerbody has a shape as seen in 
Figure 5 (for parabolic dependence). If only r_(x) sat- 
isfies these conditions, and r + (x c ) is not zero, a Screw- 
driver centerbody is formed. Note that these piecew-ise 
smooth contours are limited cases of the smooth con- 
tours described in the section 1.1.3 if power coefficients 
p and q in equations (1.4, 1.5) limit to infinity. 

For such a configuration two expanded flows near 
the nozzle wall flow into two neighboring cavities to 
meet each other at some angle a, mutually penetrate 
and more effectively mix. A flow impulse on the lat- 
eral area of the convexities increases the resulting nozzle 
thrust. 

1.5 2D Screwdriver Shaped Surfaces. We 

expect to obtain similar effects in the case of a 2D nozzle 
through the application of linear surfaces. Such a nozzle 
is a 2D analogue of an axisymmetric Screwdriver cen- 
terbody. An individual element of such a construction 
is shown in Figure 6a. The nozzle contains a rectangu- 
lar throat. A 2D plug is made with straight lines, which 
join points on the throat horizontal side with the verti- 
cal interval at the nozzle (or plug) edge. Such elements 
can be multiply repeated and form a cellular nozzle. An 
example of a two-cellular nozzle is show n in Figure 6b. 
As in the previous case of an axisymmetric flow, two 
anti-symmetric swirls form at the lower straight lines 
of two neighboring ’’ravines”. These ravines meet each 
other at the edge on the vertical interval. Two crossing 
swirl flows create large vorticity downstream. Schemat- 
ically, such swirl flows are illustrated in Figure 6c. The 
surface geometry and the surface equations are shown 
in Figure 6d. 

The proposals described in this paragraph are only 
preliminary observations, based on the approximate es- 
timations, on our experience, and some numerical sim- 
ulations. Determination of Screwdriver shape efficiency- 
can be drawn only after detailed theoretical research, 
numerical simulation, and experimental tests. 

II. Acoustic Tests of Screwdriver Shaped and 
Axisymmetric Centerbodies (Plugs) 

For examination of the described above concept , 
several experimental and numerical tests w-ere conducted. 
Acoustic tests were made in the anechoic chamber of t he 
Central AeroHydrodynamics Institute (TsAGI. Moscow ) 
[5 6] The interior dimensions of the facility within 
the wedged tips are 96x5.3x4. 2m high. Three noz- 
zle designs with exhaust jets were tested. All designs 
had axisymmetric external nozzles with a design ex- 
haust Mach number M, = 1.5 and three different internal 


parts: a) the coaxial 4-petal Screwdriver shaped cen- 
terbody/plug (SCR), b) (he coaxial axisymmetric plug 
with the length and cross sections equal to the Screw- 
driver shaped plug (CON), and c) the reference round 
nozzle (REF) without a plug with a mass flow rate equal 
to that of the previous co-annular nozzles. The two pho- 
tos of the first designs, a) and b), are shown in Fig. 7, 
and the meridional plane view with the main sizes for 
the CON design are in Fig. 8. 

The nozzle pressure ratio (NPR) is varied in the 
interval 7r c =1.6-4.5 with a cold, supersonic, and un- 
derexpended exhaust jet. The measurement procedure 
in the anechoic chamber and the method of automatic 
data processing are illustrated in Fig.9. Microphones 
(model 4136/BrulftK’err Co) with cathodic recapitu- 
lators (model 4633) are established on the meridonal 
plane at the circular arc with the radius R m = 2m 
with different observation angles 6 to a positive (down- 
stream) jet axis into an interval 30° < 6 < 105°. Micro- 
phone signals are transmitted to the magnetic recorder 
’’Sony KS-616U” through an amplifier (model 2608). 
Decoding of acoustic pressure pulsations was conducted 
by an analyzer (model 2032/Brul’ftK’err Co.) and PC- 
286 computer which provided a narrow band spectrum 
with a band size ~ A/=32Hz. A Pentium computer 
was used for transformation of narrow band spectrums 
to 1/3-octave spectrums. 

These experiments have shown an essential acous- 
tic benefit ~10-13’/. with the application of the co-annular 
nozzles by comparison with the reference round noz- 
zle. In Fig. 10 narrow band spectra of jets issuing from 
the three supersonic nozzles are compared with the the 
same mass flow rate and with the nozzle pressure ratios 
b8, 2.0, 2.4, 3.0, and 3.5. Results for the 
round nozzle without a centerbody (REF) are shown 
by black lines, the nozzle with the axisymmetric center- 
body (CON) by red line, and the nozzle with the Screw- 
driver shaped centerbody (SCR) by blue lines. Here the 
observation angle is 0 = 30°. 

The expected acoustic benefits of the Screw’driver 
shaped plug application were not obtained by compar- 
ison with an axisymmetric plug. Moreover, for some 
angles (6 = 90°), a noise increase was observed (~1- 
2/,). In Fig.ll, narrow band spectra of jets issuing 
from the two supersonic nozzles with the same mass 
flow’ rate; nozzle pressure ratios, NPR=1.6, 2.4, 3 5 
and 4.5; and three observation angles, 6 = 30°, 60°, 
and 90° are compared. Results for the CON design are 
shown by blue lines and results for the SCR design by 
yellow hues. 1/3-octave band spectra for the previous 
nozzle pressure ratios, and the four observation angles 
0 = 30°, 45°, 60°, and 90° arc in Fig. 12. Results for the 
CON design are designated by blue lines and results for 
the SCR design by yellow lines. Influence of the az- 
imuthal angle <1- on narrow band spectra at the three 


observation angles, 0 = 30°, 60°, and 90°, the three az- 
imuthal angles, * = 0*,22.5<\ and 45°, and the nozzle 
pressure ratio, NPR=3.5, is shown in Fig. 13. This illus- 
trates an inessential dependence of noise on azimuthal 
direction. Small displacement of the SCR centerbody 
upstream, so that Screwdriver shape is started exactly 
from the throat cross section, did not reduce noise and 
produced only slight changes in the 1/3 octave band 
spectra. This is illustrated in Fig. 14 where SCR corre- 
sponds to the standard centerbody location, as is shown 
in Fig. 7a, and SCR-IN corresponds to the SCR design 
moved upstream. Comparisons are made for the two 
observation angles 0 = 30° and 60°, and two nozzle 
pressure ratios, NPR=1.6 and 3.5. 

These experimental tests will be continued with the 
goal to obtain better acoustic results. In particular, we 
hope to improve acoustic characteristics by moving the 
centerbody into the nozzle and using penetrable walls 
for the SCR and/or the main nozzle. Preliminary re- 
sults for such an approach are very promising. 

III. Acoustic Tests of Bluebell Nozzles 

^Nine different Bluebell shaped nozzles were made 
at NASA LaRC and tested earlier, including the base- 
line round nozzle. The list of these nozzles, assigned 
numbers #1 thru #9, is represented in Table 1 along with 
their geometric parameters n = n p = n c ,e,6 0 . Detailed 
information on Bluebell nozzle geometry is contained 
m several papers (see, for example, [3]). Note only that 
this geometry is defined by the three main parameters: 
n — rip — n c ,e l and 6 a . These give the frequency and 
amplitudes of the sinusoidal change of the nozzle lip-line 
edge denoted as ’’shevrons” (or petals) and the nozzle 
cross section shapes are denoted as ’’corrugations”. 


Table- 1 


No. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

Ab. 

8C 

8S 

8S 

8S 

4C 

4S 

OR’ 

8S 

OR 

n 

8 

8 

8 

8 

4 

4 

0 

8 

0 

€ 

0.23 

0.23 

0.38 

0.15 

0.45 

T).45 

0 

0.23 

0 

6 0 

0.20 

0.00 

0.00 

0.00 

0.20 

0.00 

0 

0.10 

0 


In Table 1, nozzle #7 is the round convergent base- 
line nozzle and nozzle #9 is the round convergent di- 
vergent baseline nozzle. Five of these nozzles, #1,2, 5, 6, 
and 9, were sent to TsAGI, Moscow. All of the noz- 
zles have a throat radius r.=0.G37in (3.62cm) and the 
round nozzle has the exit radius r"=0.694in (1.70cm). 
The Bluebell nozzles have exit radii defined by equa- 
tion (1.2 1). The length of the convergent pari of all 
nozzles is *.=2.05in (5.2cm) and the entire length of 
the baseline round nozzle is x;=4.25in (10.8cm). All 
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of the Bluebell nozzle deisgns are based on the second 
formulation, with the use of additional sheets. ieir 
lengths are varied from the baseline nozzle by an ad- 
ditional 1.0 to 1.1 inches depending on the geometrical 
parameters selected for the nozzle. 

Abbreviations are introduced in the table for the 
tested nozzles: 8C and 8S are the 8-petal Blubell nozzles 
with corrugations (6. > 0) and without corrugations 
(’’skew” , 6 0 = 0) respectively, and 4C and 4S are th 
same for the 4-petal Bluebell nozzles. Note, that dur- 
ing the tests at NASA SAJF, in the supersonic regime 
the nozzles were operated slightly underexpanded at the 
nozzle pressure ratio NPR=4-0 and jet total tempera- 
ture T = 350°F. This corresponds to a fully expanded 
jet Mach number of 1.56 with corresponding exit ve- 
locity of 1740 ft. /sec (531m/sec). All nozzles, except 
convergent nozzle 7, were designed for a fully expanded 
Mach number M e = 1.5. This underexpanded condition 
was selected because it was expected that the nozzles 
with corrugations would reduce noise more effectively at 
this condition. In the subsonic regime the nozzles were 
operated at the nozzle pressure ratio NPR-127 and jet 
total temperature of 538"F. This corresponds to an ex- 
haust Mach number M e = 0.6 with corresponding exit 
velocity 900 ft./sec (274m/sec). At either supersonic or 
subsonic nozzle operating conditions, the nozzles were 
sized to operate at the same ideal thrust. 

The TsAGI acoustic tests were conducted with con- 
secutive increases in total pressure, in the interval 1.6-S-3-5, 
that changed the jet issuing regime. In the interval 
of nozzle pressure ratios, 1.0<NPR<1.89, a subsonic 
regime is always realized. Recall, that the reference 
round nozzle #9 is a convergent-divergent nozzle de- 
signed for a fully expanded exit Mach number, M e — 
1.5, with the rated nozzle pressure ratio, NPR—3.6 . 
Therefore, in the overcritical nozzle pressure ratios in- 
terval 1 89<NPR<3.67, the gas flow is supersonic down- 
stream of the nozzle throat. Static pressure increase 
then produces some oblique shock waves with a regu ar 
reflection from the nozzle-jet axis or with an irregular 
reflection forming a Mach disk at the axis. Shock waves 
are formed outside the nozzle creating a jet barrel-shaped 
structure. Stabilization of the shock waves inside a di- 
vergent nozzle portion is usually not observed. In an 
interval of NPR>3.67, the jet is underexpanded and 
the barrel-shaped structure is present. 

The main results of the TsAGI Bluebell nozzle acous- 
tic tests are collected in Fig. 6. Narrow band spectra of 
jets issuing from the five supersonic nozzles with the 
same mass flow rate are presented with the three noz- 
zle pressure ratios: NPR = 1.6. 2.4, and 3.5, and three 
observation angles, 0 = 30°, 60". and 90". Acoustic 
data for the reference round nozzle designed for an exit 
Mach number, M, = 15, are shown by black lines and 
denoted as OR. (nozzle #9 from the table). Data for 


the 8-petal Bluebell nozzle (BN) without corrugations 
8S (#2) are denoted by blue lines ; data for the 8-petal 
BN with corrugations. 8C (#1), by yellow lines; data 
for the 4-petal BN without corrugations, 4S (#b), by 
green lines; and data for the 4-petal BN with corru- 
gations, 4C (#5). by red lines. These spectral density 
graphics are drawn using the aforementioned analyzer 
(model 2032/Brul’ftK’err Co.) and PC-286 computer 
which provided a narrow band spectrum with a ban 
size ~ A/=32Hz. The considered frequency interval 
was from 0kHz to 25kHz, although the measurement 
system allows an interval from 0kHz to 100kHz. 

In general, the TsAGI acoustic test results for cold 
nondesigned jets have been confirmed by previous NASA 
tests showing that Bluebell nozzles provide some acous- 
tic benefits, especially in the direction of maximum noise 
production (for an observation angle, 6 - 30 ). for 
the subsonic regime, NPR=1.6, these benefits are ob- 
served in all observation angles for all Bluebell nozzles, 
excluding nozzle #5 (4C) in an interval [0>5]k z * e 
best design in this regime is nozzle #2 (8S). Recall tha 
comparison in the subsonic regime is made with the ref- 
erence round nozzle, #9, designed for the supersonic exit 
velocity with Mach number, M e = 1-5. In accordance 
with quazi-one-dimensional inviscid theory, for sue a 
nozzle with nozzle pressure ratio, NPR=l-6,^the exit 
Mach number would be equal to M t = 0.8477. PP 
cation of this theory for a Bluebell nozzle is genera y 
not valid, especially near the corrugated exit. However, 
in these cases, the flow regime remains subsonic. 

In the overcritical regime for nozzle pressure ratio, 
NPR=2.4, the reference round nozzle is overexpended. 
Strong shock wave formation in the exhaust jet leads 
to the presence of narrow discrete component (acoustic 
picks) in these spectra for all three observation angles 
It is remarkable that these picks are not observed for all 
Bluebell nozzles. It is hoped that this is not the result o 
errors connected with the imperfect data processing ap- 
paratus. Similar acoustic picks were observed in other 
experiments with the round nozzles. We can conclude 
that this shock wave system was destroyed by using 
corrugated nozzle shapes. This effect occurred for both 
designs with and without corrugations. Unfortunately, 
a fully corrugated design without petals was not made 
and it can not be definitely determined which feature 
is particularly responsible for this effect, the presence 
of petals or corrugations. This will be studied further 
in the future. However, it may be asserted that the 
presence of only petals at the nozzle exit leads to the 
elimination of these picks in the narrow band spectr 
The overall pressure level acoustic benefits of the Blue- 
bell nozzles are essentially smaller in t his regime Ilian 
in the previous subsonic regime. The best designs are 
the 8 and 4-petal corrugated nozzles, #1 and #5. 

In the near-designed regime of the issuing jet willi 



nozzle pressure ratio, NPR=3.5, acoustic benefits are 
observed only for the observation angle. 0 = 30°. In this 
case, very weak shock waves are formed downstream of 
lh^ nozzle exit and the narrow band spectra do not 
show any acoustic picks which were observed in the 
previous regime. Moreover, for the observation angle, 
@ = 90°, an increase in the noise level is observed. This 
result was surprising because the analogical NASA tests 
with closed conditions for slightly underexpanded su- 
personic hot jets have shown essential acoustic benefits 
from Bluebell nozzle application in the spectra interval 
[0,20]kHz. 

Thus, analysis of the Bluebell nozzle acoustic tests 
with cold jets at TsAGI (Moscow) have shown some ad- 
ditional effects which were not observed earlier during 
NASA LaRC acoustic tests with hot near-designed jets. 
Simultaneously, these tests bring to light some issues 
which require resolution. In particular, fabrication of 
an additional small scale Bluebell nozzle with corruga- 
tions and without petals, as well as the performance of 
identical tests with cold jets at NASA LaRC and TsAGI 
(Moscow) for proper comparison of results. 

Numerical Simulation Results 
Aeroperformance effects and gas dynamics flow char- 
acteristics were analyzed numerically. The nozzle thrust 
calculation were based on a full Navier-Stokes equa- 
tions solver (NSE), and both full and marching Euler 
codes: CFL3D [7], CRAFT [8], and Krayko-God unov 
[9]. Grid preparation and optimization was conducted 
using GRIDGEN and our own codes. The main re- 
sults were obtained using 2D and 3D versions of the 
CFL3D code which allows the simulation of both in- 
viscid and viscous flows. This code is described in de- 
tail in the CFL3D User’s Manual (Version 5.0) [10]. In 
accordance with this manual introduction ([10], page 
1): ...CFL3D (Version 5.0) is a Reynolds-Averaged 

thin-layer Navier-Stokes flow solver for structure grids... 
CFL3D solves the time- dependent conservation law 
form of the Reynolds-averaged Navier-Stokes equations. 
The spatial discretization involves a semi-discrete finite- 
volume approach. 

Upwind-biasing is used for the convective and pres- 
sure terms, while central differencing is used for the 
shear stress and heat transfer terms. Time advancement 
is implicit with the ability to solve steady or unsteady 
flows. Multigrid and mesh sequencing are available for 
convergence acceleration. Numerous turbulence mod- 
els are provided, including 0-equation, 1-equation, and 
2-equation models. Multi-block topologies are possible 
with the the use 1-1 blocking, patching, overlapping, 
and embedding. CFL3D does not contain any grid gen- 
eration software. Grids must, be supplied extraneously.” 

In Figures 16-23 some numerical results are pre- 
sented. In Fig. 16a, b Cartesian (x,y,z) and cylindri- 
cal (;,r. 0) coordinate systems and a typical 3D grid 


are shown with indicies ( i,j,k ) running from (1,1,1) to 
(idim, jdim,kdim) in a normalized rectangular 3D re- 
gion. This grid was used for numerical simulation of the 
supersonic nozzle with the 4-petal Screwdriver shaped 
centerbody (plug). The region of direct calculation is 
limited by 6 boundaries with the following boundary 
conditions: 1) In the 0 direction by the two planes of 
symmetry 0 = 0° (r=0) and 6 = 45° with the bound- 
ary conditions for symmetry planes (or conditions for 
an inviscid flow on a solid surface). 2 ) In the x direc- 
tion by the two cross section planes, x=0 (inlet), and 
x = x e (exit). At the inlet so called "inflow /outflow” 
([8J) inviscid boundary conditions are used which define 
all of the non-dimensional flow parameters if pressure 
and temperature ratios are given. Here, ID Riemann 
invariants are applied in the direction normal to the 
boundary. At the exit, extrapolation boundary condi- 
tions are used which are cell-center type boundary con- 
ditions. 3 ) In the radial r direction, either inviscid or 
viscous solid wall conditions are imposed. 

Initially, the numerical calculations focused on the 
simulation of the inviscid and viscous flow issuing from 
jets having axisymmetric and Screwdriver centerbod- 
ies. However, these geometries provide several features 
which make performing such calculations difficult. There- 
fore, computations were performed which focused on 
the nozzle portion of the flow to gain some experience 
and provide some understanding of the issues involved 
in the full jet simulations. While many of the results are 
preliminary, it is useful to study some results and con- 
clusions. Figure 17 illustrates some results for viscous 
flows (using Menter’s k-u SST turbulence model [11]). 
within nozzles having axisymmetric and Screwdriver 
centerbodies (nozzle pressure ratio NPR = p 0 /Poo = 
2.296 and nozzle temperature ratio NTR = T„ IT X = 
2.817). Figures 17b and 17c show the mach contours 
for the axisymmetric centerbody and the 0 = 45° plane 
for the Screwdriver centerbody. The general profile of 
the two flows is similar, as also seen in figures 17d and 
17e, which provide some indication of the flow near the 
surface of the centerbody. 

One of the issues of interest was the introduction 
of vorticity to the nozzle flow by the Screwdriver cen- 
terbody, as opposed to a simple axisymmetric center- 
body. Figures 18a, b, and c provide some illustration of 
the general viscous flow within the nozzle. Figures 18d 
and 18e show the mach contours at the nozzle exit for 
an inviscid and viscous case. Superimposed on these 
contours are flow streamlines which show a consider- 
able change in the flow field for the inviscid and viscous 
cases. While some of this vorticity may be generated 
by the numerical scheme itself, it would be hoped that 
some vorticity is actually generated by the presence of 
the Screwdriver centerbody, leading to enhanced mix- 
ing. 


0 


la figure 19, the results from studying the influence 
of the axisvmmctric centerbody location on the inviscid 
(low field within the nozzle are presented. In looking at 
the mach contours, it may be seen that movement of 
the centerbody toward the throat produces shock waves 
which are reflected downstream. In this case, however, 
this strong shock structure does not greatly effect the 
thrust production. 

With experience from the nozzle calculations in 
hand, numerical simulations of some full jet flows were 
performed. Figure 20 shows results for an inviscid jet, 
with an axisymmetric centerbody, employing two dif- 
ferent computational grids. In this figure, grid 1 corre- 
sponds to a coarse grid with dimensions 162 x 66 in the 
i and j directions, respectively. Grid 2 is a much finer 
grid with dimension 306 x 154. The residual history 
shows that the coarse grid calculations converges much 
better than the fine grid calculation. The comparison of 
the mach contours seem to indicate that this behavior 
may be the result of problems at the outflow boundary. 

In this region, a large flow gradient is present which 
may lend to some instability. 

Simulation results for a jet with a Screwdriver cen- 
terbody are included in figure 21. It appears that there 
is more mixing in this case. However, as in figure 20, 
it may be seen that the solution is not fully converged. 
This may again be caused in part by difficulties at the 
jet outflow boundary. This case may however provide 
some insight into the form of the flow field. Figure 22 
shows mach contours and streamlines for this inviscid 
flow on the surface of the Screwdriver centerbody. It 
can be seen that this flow incorporates flow features 
which would not be present with an axisymmetric cen- 
terbody. 

One final area of interest is the effect of the cen- 
terbody geometry on thrust. Figure 23 presents the 
drag produced by the subsonic portion of the nozzle 
minus the thrust produced by the centerbody. Study- 
ing this “effective” drag provides an indication of the 
thrust production with various centerbody geometries. 
The preliminary conclusions of this numerical analysis 
indicate that the nozzle with the SCR plug provides less 
thrust than the nozzle with the impenetrable CON plug 
by a factor of -1-2’/.. Detailed graphic visualization of 
the gas flow in the nozzle and exhaust jet have revealed 
some important acoustic effects in the experiments con- 
ducted. With the experience gained from these prelim- 
inary calculations, it is hoped that further numerical 
simulations will provide insight into the study of these 
jet flows. 

Conclusions 

Several different corrugated designs, including Screw- 
driver shaped centerbodies and Bluebell nozzles were 
tested experimentally and numerically. Acoustic tests 
conducted at TsAGI, Moscow and numerical simula- 


tions at NASA LaRC have provided preliminary re- 
sults which indicate that the tested Screwdriver shaped 
centerbody did not produce the expected acoustic and 
aeroperformance benefits. Acoustic tests performed on 
the Bluebell nozzles essentially confirmed, in general, 
previous NASA LaRC experiments which show acoustic 
benefits from these designs (including subsonic regimes). 
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Fig. 1 


Fig.l A Screwdriver centerbody (plug), a) -3D lateral view: 1 - nozzle, 2 -centerbody, 3 -upstream swirls 
at the line of symmetry, 4 -downstream swirl behind a centerbody along the axis of symmetry; b) - up- view 
from downstream. 
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DIFFEREN T MODIFICATIONS OF SCREWDRIVER 

CENTERBODIES 


nozzle 





, >V\ 










Fig.2 Different modifications of screwdriver centerbod- 
ies with symmetrical subregions: a-conical centerbody; 
b- two symmetrical subregions, n=2; c-n=3; d-n=4; 
e- n=6- f-Cartesian coordinates and a surface equation. 
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Fig.3 A Sonic nozzle with the Screwdriver Centerbody as a unit. The model was drafted for 
experimental tests at the TsAGI, Moscow. 

a) common view: 1-sonic nozzle throat, 2-Screwdriver centerbody, 3-external nozzle surface, 4- plate to keep 
the holders, 5-bolts to keep the plate and to compress the holders, 6-plane ground for the plate; 

b) view without the external cone; c) lateral view without the external cylindric part; d) back view; e) internal 

view without some shades; f) axisymmetric centerbody with the equal cross section areas as well as the 
Screwdriver centerbody. ia 






Fig. 4 


Fig.4 4-petal Screwdriver centerbody (plug). 1 - nozzle, 2 -centerbody. a) -Srewdriver centerbody with the 
arc-shaped edge, symmetrical location of petals without a lift and momentum, b) -Screwdriver centerbody 
(plug) with a rotative momentum. 
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Chisel centerbody (afterbody) geometry. 



n c =4 


Fig • 5 


Fig. 5 A Chisel centerbody (plug) geometry. 1) n c = 8: a)-3D view, 1- nozzle, 2- centerbody. b)- centerbody 
cross section. 2) n c = 4: c)-3D view, 1- nozzle, 2- centerbody. d)- centerbody cross section. 
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Fig.6 A 2D analogue of an axisymmetric Screwdriver centerbody (or a plug), 
a) -a single 2D element of the nozzle (n^l). A surface of the supersonic part is linear: the yellow element 
contains straight lines (red), which join points of the throat horizontal side with the vertical interval at the 
nozzle (or the plug) edge; the two white elements contain straight lines (blue), which join points of the throat 
vertical sides with the two horizontal half-intervals at the nozzle (or the plug) edge, b) - a two-contour element 
(n e =2). c)-Schematic swirls into a 2D Screwdriver nozzle cell flow; d)- surface geometry. 
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Fig. 7 Axisymmetric convergent-divergent nozzle with the Screwdriver-shaped (a) and ax- 
isymmetric (b) centerbodies. Both centerbodies have the same areas at the cross section 
located on the same distance, x = £j=const, from the throat, x = 0. These designs were 
tested in the anechoic chamber at the TsAGI, Moscow. 

14 




Fig.8 The picture of the nozzles with the Screwdriver-shaped (SCR) and axisymmet- 
ric (CON) centerbodies (a). The main sizes (mm.) of the nozzle CON (b). 



Microphone B&K type 4135 and 4136 
•Preamplifier BK type 2633 

•Amplifier B&K type 2608 
assette recorder Sony KS-616U 


^Signal analyser B&K type 2032 


— -Computer * 


Narrowbands spectra 


1/3 -octave spectra 


Fig. 9 The scheme measurement in the anechoic chamber AK-2 at the 
TsAGI, Moscow, and the scheme of automatic data processing. 
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Fig-11 Narrowband spectra of jets issuing from the two supersonic nozzles with the same 
mass flow rate, with the four nozzle pressure ratios: NPR= 1.6, 2.4, 3.5 and 4.5, 
and three observation angles, 0 = 30°, 60°, and 90°. Experimental results for the 
round nozzle with the axisymmetric centerbody (CON) lire shown by blue lines, for 
the nozzle with the Screwdriver-shaped centerbody (SCR) are shown by yellow lines. 
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Fig. 12 1 /3-octave band spectra of jets issuing from the two supersonic nozzles with the 

same mass flow rate, with the four nozzle pressure ratios: NPR= 1.6, 2.4, 3.5 and 
4.5, and four observation angles, 0 = 30°, 45°, 60°, and 90°. Experimental results 
for the round nozzle with the axisymmetric centerbody (CON) are shown by blue 
lines, and for the nozzle with the Screwdriver-shaped centerbody (SCR) are shown 
by yellow lines. 
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Fig 13 Influence of the azimuthal angle S on narrow- Fig.14 Influence of location of the Screwdriver-shaped 

band spectra at the three observation angles, 0 = 30°, 00° centerbody on 1/3-octave band spectra for the two ob- 

and 90°, and nozzle pressure ratio, NPR=3.5. servation angles, © = 30° and 60°, and two nozzle 

pressure ratios, NPR=1.6 and 3.5. SCR-corresponds to 
standard location as is shown in Figure 7a, and SCR- 
19 IN-corresponda to the Screwdriver-shaped centerbody 
moved upstream so that Screwdriver shape is started 











BLUEBELL NOZZLES ACOUSTIC TESTS: (TsAGI, MOSCOW, CHAMBER AK-2) 
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3-D grid for jet flow simulation issuing from the nozzle 
with the screwdriver centerbody 


a) 3-D grid view 


y b) 2-D arid at x=0.0 



Y c) 2-D grid at x=8.5 



2.0 1.5 1.0 0.5 0.0 



Fig. 16 3-D grid for jet flow simulation issuing from the nozzle with the screwdriver centerbody. a) 3-D grid 
view: in the XOY-plane, a plane of symmetry, 0 = 45°, and on the centerbody’s surface, b) 2-D grid at the inlet, 
x=0.0, where inlet flow conditions are specified, c) 2-D grid downstream of the end of the centerbody, x =8.5. 
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a) Residual history 



b) Mach contours: axisymmetric centerbody 
2.0 r (Re = 1.E6) 
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c) Mach contours: screwdriver centerbody 
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d) Flow near axisymmetric centerbody 



e) Flow near screwdriver centerbody 
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Fig. 17 Residual history and Mach contours for nozzles with axisymmetric and screwdriver centerbodies (nozzle 
pressure ratio NPR = p 0 /Poo = 2.296, nozzle temperature ratio NTR = T 0 /Too = 2.817). a) Residual history 
for the two cases, b) and c) Results of Navier-Stokes simulations with the two types of centerbody. d) and e) 
Illustration of the flow profile near the centerbody surfaces. 
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Mach Contours for viscous and inviscid nozzle flow simulations 



Fig. 18 Mach contours for viscous and inviscid flow simulations (nozzle pressure ratio NPR = po/poo = 2.960, 
nozzle temperature ratio NTR = T 0 /Too = 2.817). a) 3-D view of mach contours for a viscous flow within a nozzle 
with a screwdriver centerbody. b) Mach contours in the 9 = 0° plane ( z = 0). c) Mach contours in the 9 = 45° 
plane, d) Mach contours with superimposed streamlines at the nozzle exit (end of the screwdriver centerbody) 
for an inviscid flow case, e) Mach contours with superimposed streamlines at the nozzle exit for a viscous flow 
case. 
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Fig. 19 Influence of axisymmetric centerbody location on the inviscid flowfield within a M=1.5 designed nozzle 
(nozzle pressure ratio NPR = p 0 /Poo = 2.96). a) Mach contours for centerbody located near x — 2.0. b) Mach 
contours for centerbody located near x = 1.0 (closer to the nozzle inflow). 
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a) Residual history 



b) Mach contours: Grid 1 



Fig* 20 Residual history and mach contours for an inviscid jet with an axisymmetric centerbody employing 
two different computational grids (nozzle pressure ratio NPR = p 0 /p<x> = 3.67). a) Residual history, b) Mach 
contours for grid 1 (coarse grid with dimensions: 162 x 66). c) Mach contours for grid 2 (fine grid with dimensions: 
306 x 154). 
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a) Residual history 
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b) Mach Contours in the 0 = 0° plane 
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c) Mach Contours in the 6 = 45° plane 
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Fig. 21 Residual history and mach contours for an inviscid jet with a screwdriver centerbody (nozzle pressure 
ratio NPR = p 0 /Poo = 2.296). a) Residual history, b) Mach contours in the 0 = 0° plane (z = 0.0). c) Mach 
contours in the 0 = 45° plane. 
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Fig. 22 Mach contours and streamlines on the surface of the screwdriver centerbody of an inviscid jet flow 
(nozzle pressure ratio NPR = p 0 /poo = 2.296). 












Drag history comparison 
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Fig. 23 Drag history for inviscid and viscous nozzle flows with axisymmetric and screwdriver centerbodies 
(nozzle pressure ratio NPR = p c /p x = 2.960, nozzle temperature ratio NTR = T 0 /Toa = 2.817). a) Comparison 
of several cases showing progression toward a converged value, b) Comparison of of viscous nozzle flows with a 
screwdriver centerbody employing grids with different stretching near the centerbody surface. (Grid 1: grid used 

in Euler calculations, nearly uniform; Grid 2: grid employing stretching near the wall in attempt to capture some 
boundary layer effects). 
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